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We study the beam spin asymmetry A 



sin <j)f, 
LU 



in semi-inclusive tv electroproduction contributed 



by the T-odd twist-3 distribution function g ± (x,k^). We calculate this transverse momentum 
dependent distribution function for the u and d quarks inside the proton in a spectator model 
including the scalar and the axial- vector diquark components. Using the model results, we estimate 
the asymmetry A s ^ in the ep — > e'-K°X process in which the lepton beam is longitudinally 
polarized. The model prediction is compared with the data measured by the CLAS and HERMES 
Collaborations, and it is found that our numerical results agree with experimental data reasonably. 
Especially, our results can well describe the CLAS data at the region where the Bjorken x and the 
pion transverse momentum is not large. We also make a prediction on the asymmetry A a ^/' h in ty° 
electroproduction at CLAS 12 using the same model calculation. 



PACS numbers: 12.39.-x, 13.60.-r, 13.88.+e 



I. INTRODUCTION 



Understanding the origins of the single spin asym- 
metries (SSAs) appearing in high-energy semi-inclusive 
processes is one of the important goals of QCD spin 
physics Substantial SSAs have been measured by 

the HERMES Collaboration g|-G3 , the COMPASS Col- 
laboration [iMl, and the Jefferson Lab (JLab) [l4Hl8| 
in semi-inclusive deep inelastic scattering(SIDIS). It is 
found that the T-odd transverse momentum dependent 
(TMD) parton distribution functions (DFs) (Ta^ or 
fragmentation functions (FFs) [25[, under the TMD fac- 
torization j2o| framework, play central roles in the ob- 
served SSAs. The T-odd TMDs describe the correlations 
between the transverse motion of the parton and its own 
spin or the spin of the initial-state hadron, thereby encod- 
ing much richer information about the partonic structure 
as well as the QCD dynamics inside hadrons than what 
can be learned from the collinear DFs. 

In the particular case of a longitudinally polarized 
beam colliding on an unpolarized target, an asymme- 
try with sin<ph modulation, the so-called beam SSA, 
emerges. The CLAS Collaboration [3 Q3, 01 at JLab 
and the HERMES Collaboration Q have measured this 
asymmetry in pion electroproduction in the magnitude of 
several percents that cannot be explained by perturba- 
tive QCD 27]. Different mechanisms have been proposed 
to generate such asymmetry, such as the Boer-Mulders 
effect [28[ and the Collins effect [29|, |30(, involving the 
chiral-odd distribution or fragmentation functions. In 
Refs. p3lH33j . a new source contributing to the beam 
SSA has been identified, either from the model calcu- 
lations [3l], HII or from the updated decomposition of 
the unpolarized quark-quark correlator [33| . where the 
twist-3 TMD g ± (x,k^,) plays a crucial role. As a T- 



odd chiral-even TMD, can be regarded as an analog 
of the Sivers function [l9j at the twist-3 level, because 
both of them require quark transverse motion as well as 
initial- or final-state interactions via soft-gluon 

exchanges to receive nonzero contributions. Therefore, 
studying beam SSAs may provide a unique opportunity 
to unravel the role of quark spin-orbit correlation at twist 
3. 



In this work, we present an analysis on the beam SSA 
in neutral pion production, based on the effect 
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from g- 1 . We calculate the function g ± (x,k^) for the 
u and d quarks inside the proton, using the spectator 
model with scalar and axial-vector diquarks. Different 
types of spectator model have been widely used to calcu- 
late TMDs for the nucleon [H |34H4^ and the pionjil- 
[46| . We will adopt a specific model given in Ref. |4l|. in 
which the authors consider the isospin of vector diquarks 
to distinguish the isoscalar (ud-like) and isovector (mi- 
like) spectators. Furthermore, in that model the free pa- 
rameters are fixed by reproducing the parametrization of 
unpolarized and longitudinally polarized parton distribu- 
tions. The above-mentioned feature of the model allows 
us to perform the phenomenological analysis in a deep 
sense. Using the calculated g ±u and g ±d , we estimate 
the beam SSA A s L "fj }h in neutral pion production at the 
kinematics of CLAS and compare our results with the 
CLAS data [H measured recently with high precision. 
We also make a comparison between our calculation and 
the 7T° data measured by the HERMES Collaboration Q 
for further testing. Finally, we present the prediction of 
Af^ 4 " 1 for 7T° production at CLAS12. 

The paper is organized as follows. In Sec. [Hi we present 
the details on the calculation of g- 1 in the spectator model 
with an axial-vector diquark, and discuss the flavor, x 
and kr dependencies of g , In Sec. IIII1 we analyze the 
beam SSA A L "^' h in ir° production numerically at CLAS, 
HERMES, and CLAS12 based on the calculated g 1 - for 
the u and d quarks. We address our conclusion in Sec. lIVI 
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II. p x OF THE PROTON IN THE SPECTATOR 
MODEL WITH AN AXIAL- VECTOR DIQUARK 

In this section, we present the detailed calculation on 
the T-odd twist-3 TMD g ± (x, kj-) of the proton in a spec- 
tator model with axial- vector diquark. The starting point 
of the calculation is the gauge- invariant quark-quark cor- 
relator for the unpolarized nucleon 

$W(z,fc T )= / %- ^ e' fe - g (Pfe-(0)£[Q-,oo-] 



.., (2tt)3 " ^ ,TJ ' 
xC[0 T ^ T }C[^-,C}UO\P), (i) 

where [+] denotes that the gauge link appearing in $ 
is future-pointing, corresponding to the SIDIS process. 
As an interaction-independent twist-3 distribution, g 1 - 
naturally appears in the decomposition of the quark- 
quark correlator at the subleading order of 1/P + ex- 
pansion, after including the light-cone vector n_ that 
defines the direction along which the path-order expo- 
nential runs [H E3 



$ [+1 (x,fc T ) 



= J &-$ [+1 (F,fc;n_) 

_ M ( ± e P Tl p k To 
~ 2P+\ 9 75 M 



(2) 

here, • • • stands for the other twist-3 TMDs that are not 
taken into account in this paper. The distribution g , 
therefore, can be deduced from $[ + l(a;, kr) by taking the 
trace with proper Dirac matrices 



P 



^(*,4) = Tr[$M 



7 75j 



(3) 



Since g is a T-odd DF, we need to consider the effect 
of the gauge link to generate a nonzero result, in analogy 
to the Sivers function and the Boer-Mulders function. 
However, calculations on g^ based on the scalar diquark 
model (with pointlikc proton-quark-diquark coupling) as 
well as the quark-target model [48[ show that g ± and the 
other twist-3 T-odd TMDs suffer light-cone divergence. 
This feature is in contrast to the case of leading-twist 
T-odd TMDs and has been recognized as a theoretical 
challenge in deriving a TMD factorization proof in SIDIS 
at twist 3 HHim. In Refs. (11,113, the authors pointed 
out that light-cone divergence can be avoided by means 
of a phenomenological approach in which form factors 
are applied for the proton-quark-diquark coupling. In 
Ref. [51(, all sixteen twist-3 T-odd TMDs were calcu- 
lated in the scalar diquark model by adopting the dipo- 
lar form factor for the proton-quark-diquark coupling. 
In this paper, we extend the calculations in Refs. [48l. l5l| 
using the spectator model with an axial-vector diquark 
to obtain g 1 - for both the u and d quarks. The model 
we apply in the calculation is the version developed in 
Ref. [4l| , which was originally proposed in Ref. [34 1 . The 
same model has been adopted to calculate various twist-3 
quark-gluon-quark correlators in Ref. 52] . The common 
feature of spectator models is that the proton with mass 
M is supposed to be constituted by a quark with mass 
to and a diquark with mass Mx , and the diquark X can 
be either a scalar one (denoted by s) or an axial-vector 
one (denoted by v). 



We perform the calculation in the Feynman gauge 
and expand the gauge link to the first order (one-gluon 
exchange) to obtain the quark-quark correlators con- 
tributed by the scalar diquark and the axial-vector di- 
quark components 



®sij(x, k T ) = -ie q 



d 4 q 

(2tt) 4 J (2^)4 



dk' 



U(P,S)T s (k 



2 (^ + m) 



k 2 — to 2 + ie 



(P - k) 2 - M 2 + ie 



r+ 



q 2 + ie (P - k + q) 2 - M 2 + ie 



+ to) 



$ V ij(x, k T ) = -ie. 



1 



d 4 q 



dk- 



9 (2tt) 4 J (2tt) 4 J q-r+ie 
d0a (P-k + q) 



(k — q) 2 — m 2 + ie 
U(P,S)TP(k 2 ) 



T s ((k-q) 2 )U(P,S) 



+ H.c, 



Qt + to) 



k 2 — to 2 + ie 



k+=xP+ 

d pa {P-k) 



(P - k) 2 - M 2 + ie 



q 2 +ie(P- k + q) 2 - M 2 + ie 



(ft — i + m) 



(k — q) 2 — m 2 + ie 
I 



-T^k-q) 2 )U(P,S) 



B.C. 



k+=xP+ 



(4) 



(5) 



where e q is the charge of the quark, T s / v denotes the and or T^ a ^ is the vertex between the gluon and the 
nucleon-quark-diquark vertex with the form [lij 



T s (fc 2 ) = . 9s (fc 2 ), T^fc 2 ) = ^|V 7 5 , (6) 
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scalar diquark or the axial-vector diquark 

T' 1 =ie s (2P-2k + q)", 
T ^ a /3 = _ i6v [(2P -2k + qYg a P - (P - k + q) a g 



(7) 



I'd 



(8) 



with e s / v denoting the charge of the scalar/axial- vector 
diquark. In Eq.([S]) we use to denote the summa- 
tion over the polarizations of the axial- vector diquark for 
which we choose the following form (53j 



d^ v {P - k) = - g^ v + 



(P - k^n-u + (P - k) v n- M 



(P -k)- n- 



Ml 



[(P-k) -n. 



2 ,( '—fi"'—i'- 



(9) 



The advantage of the above choice has been argued in 
Ref. [4l[. To obtain finite results for g^, we also choose 



the dipolar form factor for gx{k 2 ) 

k 2 — m 2 



gx(k 2 )=N x 



= N 



x- 



\k 2 -A 2 x \ 2 
{k 2 - m 2 )(l - x) 2 
{k T + L 2 X ) 2 



(10) 



for X = s, v, where Ax is the cutoff parameter, Nx is the 
coupling constant (which also serves as the normalization 
constant), and L\ has the form 

L 2 X = (1 — x)A x + xM 2 x - x(l - x)M 2 . (11) 

Performing the integrations in Eqs. Q and ([5]) over 
k~ , q + and q~ , the quark-quark correlators for the un- 
polarized nucleon in the spectator model are simplified 

as 



$ s (x, k T ) = -ie q e s N, 



2 (1-*) 3 



1 



d 2 q T + m){lf + M){ k 1 + m)\ 



$ v (x,k T ) = -ie q N? (1 X) 



32tt 3 P+ (L 2 + k 2 T ) 2 J (2tt) 2 q\ (L 2 + (k T - q T ) 2 ) 2 
2 1 f d 2 q T 



(12) 



128^ 3 (P+) 2 (L 2 +fc 2 ) 2 J (2tt) 2 



k + =xP + 

d pa (P - k) (-iT+^) d a p (P-k + q) 



+ m)l"{ if - M)Y{ # + m)} 



q 2 T {L 2 v + (k T ~q T ) 2 f 



q + =0 
k + =xP+ 



(13) 



Using Eq. (j3|) and performing the integration over gr, we 
arrive at the expressions for the distribution g from the 

1 There is a typo in Ref. [H, the factor (1 - x)A 2 + (1 + x)M 2 - 
(1 — x)(l — 2x)M 2 in the numerators of the right-hand side of 



scalar and the axial- vector diquark components 1 

Eqs. (32) and (36) should be (1 - x) A 2 + (1 + x) M 2 - (1 - x)M 2 . 



g ±s (x,k 2 T ) = - 



g ±v (x,k 2 T ) 



N 2 (l~x) 2 e s e q 
(32tt 3 ) 4tt 
2 . _ N 2 (l - x) 2 e v e q 



(1 - ,e)A 2 + (1 + x)M 2 - (1 - x)M 2 



L 2 (L 2 + fc 2 ) 3 
(1 - x){xM + m) 2 + (1 - x) 2 M 2 - Ml + xL 2 v 



(32tt 3 ) 4tt 
x 

(1 — x)/c 2 n(7v 2 -|- fc^i) 2 



(l-x)L 2 v (L 2 v + k 2 T f 



In 



L ■ 



L 2 



(14) 



(15) 



The function g for the u and d quarks can be con- 
structe d by g and g . Here, we follow the approach 
in Ref. [41j in which the two isospin states of the vector 
diquark are distinguished; that is, the vector isoscalar 
diquark v(ud) is denoted by a with mass M a and cutoff 
parameter A a , while the vector isovector diquark v(uu) is 
denoted by a' with different mass M a i and cutoff param- 



eter A a / . Hence, by applying the relation between quark 
flavors and diquark types, we can obtain the distributions 
g ±u and g ±d by the following form |4lj : 

9 ±u = c 2 g^ + c 2 a g ±a , (16) 
9 Xd = c' 2 g Xa ', (17) 
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Diquark Afx(GeV) Ax (GeV) c x Nx 

Scalar s (ud) 0.822 0.609 0.847 11.400 

Axial-vector a (ud) 1.492 0.716 1.061 28.277 

Axial-vector a' (uu) 0.890 0.376 0.880 4.091 

TABLE I: Values for the free parameters of the model 
taken from Ref. [U, which are fixed by reproducing the 
parametrization of unpolarized 54] and longitudinally polar- 
ized [5^] parton distributions. 

where c s , c a , and c a < represent different couplings that 
are the free parameters of the model. The same form 
also holds for the other TMDs. 

For the parameters Mx, Ax and cx (X — s,a, a') 
needed in the calculation, we also adopt the values from 
Ref. [4l| , as shown in the first three columns of Table. |U 
The fourth column shows the values for the correspond- 
ing normalization constants Nx {X = s, a, a'), which are 
obtained from the normalization condition for the unpo- 
larized TMD. The quark mass is chosen as m = 0.3 GeV. 
To convert our calculation to real QCD, we use the fol- 
lowing replacement for the combination of the charges of 
the quark q and the spectator diquark X 

^ -> ~C F a s (18) 

47T 

and choose a s ~ 0.3 in our calculation. The minus sign 
in the above equation comes from the fact that in QCD, 
a hadron is color-neutral; thus, the color charges e q and 
ex should have the opposite signs. 

In the left panel of Fig. [TJ we plot the functions xg^ u 
(solid line) and xg ±d (dashed line) vs kr at x — 0.3, 
while in the right panel of the same figure, we display 
the x dependence of xg ±u and xg ±d at fcy = 0.3 GcV. 
The plots in Fig. [1] correspond to the results in the scale 
= 0.3 GeV 2 , which is the scale used in Ref. [4l| to fit 
the parametrization of fi(x) [54| . As we can see from 
Fig. [D the dominance of u quark contribution is evident 
in the adopted spectator model; that is, g ±u is several 
times larger than g^~ d by size. Our results show that g ±u 
is positive for all x and kr regions; while g ±d is negative 
in the small x region and turns to be positive in the re- 
gion x > 0.15, i.e., there is a node in the x dependence of 
g ±d . Also, the kr dependencies of the u and d quark dis- 
tributions are different since g ±d approaches zero faster 
than g ±u when kx increases. 

III. NUMERICAL RESULTS FOR BEAM SPIN 
ASYMMETRY 

In this section, we will use our model resulting g 1 - 
to calculate the beam SSA A s £y h in ir° electroproduc- 
tion, as precise measurements on A s ^j' h of a neutral pion 
for different x and Pt bins have been performed by the 
CLAS Collaboration pj} at JLab recently, in SIDIS by 
a 5.776 GeV longitudinally polarized electron beam off 
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FIG. 1: Left panel: model results for xg ±u (solid line) and 
xg ±d (dashed line) as functions of kr at x = 0.3; right panel: 
model results for xg ±u (solid line) and xg ±d (dashed line) as 
functions of x at kr — 0.3 GeV . 

an unpolarized hydrogen target. Earlier, the HERMES 
Collaboration measured the beam SSAs for neutral and 
charged pions using a 27.6 GeV beam. We compare the 
results for ^ i ^f' h with the neutral pion data from CLAS 
and HERMES to test our model calculation. Then, we 
will make new prediction on the asymmetry A s ^j for 
7T° electroproduction at CLAS12 using the same model 
calculation to study the prospects to access g 1 - at CLAS 
after the 12 GeV upgrade is realized. 

The semi-inclusive leptoproduction process that we 
study can be expressed as 

e(£) + p(P) -> e'(£') + h(P h ) + X(P X ) , (19) 

where £ and £' are the four-momentum of the incom- 
ing and scattered lepton, and P and Ph are the four- 
momentum of the target nucleon and the detected final- 
state hadron h, respectively. 

The variables to express the SIDIS cross section are 
defined as 

_ Q 2 _ P-q _ P-P h _ 2Mx 

X ~2P~q' V ~P~V Z ~"P^T' 7__ Q - ' 
Q 2 = -q 2 , s = (P + £) 2 , W 2 = {P + q) 2 , (20) 

where q = £ — £' is the four-momentum of the virtual 
photon and W is the invariant mass of the hadronic final 
state. 

The reference frame we adopt here is that the virtual 
photon and the target proton are collinear and along 
the z axis, with the photon moving toward the target 
in the positive z direction, as shown in Fig. [5] We use 
kx to denote the intrinsic transverse momentum of the 
quark inside the proton for the DFs, with Pt to denote 
the transverse momentum of the detected hadron. The 
transverse momentum of the hadron h with respect to 
the direction of the fragmenting quark is denoted by p T , 
which appears in the TMD FFs. The azimuthal angle 
between the lepton and the hadron planes is defined as 
(f>h, following the Trento convention |5q |. 

Up to subleading order of 1/Q, the differential cross 
section of SIDIS for a longitudinally polarized beam (with 
helicity A e ) off an unpolarized hadron has the following 
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FIG. 2: The kinematical configuration for the SIDIS process. 
The initial and scattered leptonic momenta define the lepton 
plane (x — z plane), while the detected hadron momentum 
together with the z axis identify the hadron production plane. 



general expression [57} : 

do 2na 2 



dxdy dzhdP^,dcj)h xyQ 2 2(1 — e) 



{Fi 



+\ a y/2e{l-e)am4>h FT 1 ** 



11 

2x 



uu 



}■ 

(21) 



where Fujj is the helicity-averaged structure function, 
while is helicity-dependent structure function re- 

sulting from the antisymmetric part of the unpolarized 
hadronic tensor. The first and second subscript of the 
above structure functions indicate the polarization of 
beam and target, respectively. It is F^^ h that gives 
rise to the sin fa beam SSA. The ratio of the longitudi- 
nal and transverse photon flux is given by 



l-y + y 2 /2 + -f 2 y 2 /4' 



(22) 



In the parton model, the two structure functions in 
Eq. (l2~Tj) can be expressed as the convolution of TMD DFs 
and FFs. based on the tree-level factorization adopted in 
Ref. [13 ■ With the help of the notation 

C[wfD] =x^2e 2 f d 2 k T f d 2 p T S 2 {zk T -P T +p T ) 

xw(k T ,p T )fi(x,k 2 )Di(z,p 2 T ), (23) 

Fjjv an d Ff^^ h can be written as [13] 

F uu =C[f 1 D 1 ], (24) 



r lu ~ n 



P T k T M h U± E . 

hf h x q D\ 

M \ M z 1 



Pt Pt ( M h G± , 

11 h X eiif" 

M h \ M z 1 



(25) 



where P T = with P T = \P T \, M and M h are the 
nucleon and hadron masses, respectively. The functions 
G x and E are the interaction-dependent twist-3 FFs that 
come from the quark-gluon-quark correlator for FFs. The 
former one is T-odd, while the later one is T-even. They 
can be connected to the interaction-independent twist-3 
FFs G x @ and £ H by the following relations [57|: 



G-L _G ± m , 
~~~~ Mft 1 ' 



E E m 

— = -n-F>i- (26 

z z M h K > 



The beam-spin asymmetry A s ^j in single-pion pro- 
duction off an unpolarized target thus is expressed as 



fdxfdyfdz . x ( 1 + ^-)J2e{l-e) Ff? T + 

Jdxfdyfdz X y Qi2( £ e) x [1 + ^jFuu 



for the PT-dependent asymmetry. The x-dependent and 
the z-dependent asymmetries can be defined in a similar 
way. 

The structure function F^y receives various contri- 
butions from the convolution of the twist-2 and twist-3 
TMD DFs and FFs, as shown in Eq. ([25]). The h{E 
term, instead of the h^E in Eq. (|2"5j) . has been calcu- 
lated in Ref. (28|. The contribution from the g x D\ term 
has been studied in Refs. [3lM33l . [59l |. These two terms 
generate the asymmetry through the effects of the T-odd 
distribution functions, namely, the twist-2 Boer-Mulders 
function 



21( and the twist-3 g ± . Each distribution 



represents a specific spin-orbit correlation of the initial 
quark inside the nucleon. The beam SSA of the 7r + mc- 



effect eH^, where is the T-odd Collins FF [H 
e is the chiral-odd twist-3 DF I58ll60ll. 



and 



son has also been analyzed 29, 30] based on the Collins 



In the following, we will calculate the beam SSA in 
semi-inclusive pion electroproduction contributed by the 
g Di term. To do this, we first neglect the quark-gluon- 
quark correlators (often referred to as the Wandzura- 
Wilczek approximation [611 ]) for FFs, which is equiva- 
lent to setting all the functions with a tilde to zero. It 
is worthwhile to point out that a calculation from spec- 
tator model [13] as well as a mo del- independent anal- 
ysis [62j of the T-odd collinear quark-gluon-quark cor- 
relators shows that the gluonic (partonic) pole contri- 
butions for FFs vanish. The FF G (x,jtj,) appears in 
the decomposition of the T-odd part of the TMD quark- 
gluon-quark correlator [57l l63j , for which the gluonic pole 
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contribution should play an essential role. Whether the 
vanishing gluonic pole matrix elements for collinear FFs 
can be generalized to the case of TMD FFs deserves fur- 
ther study 64]. Nevertheless, we ignore the G 1 - and E 
contributions based on the Wandzura-Wilczek approxi- 
mation. Therefore, there are only two terms inside the 
square brackets in the right-hand side of Eq. (I25[) re- 
maining. Moreover, we consider merely the beam SSA 
of the 7T° production, since the fascinating fact that the 
favored and the unfavored Collins functions have similar 
sizes but opposite signs [H, [66[ suggests that the eH^ 
term leads vanishing beam SSA in ir° electroproduction. 
Specifically, the isospin symmetry determines that the 
7T° FF should be the average of ir + and tt~ FFs, so that 



H 



{Hi 



±fav 



Hl Unf )/2: 



0. Thus, in the follow- 



ing calculation, we can just take into account the term 
g^Di and obtain 



r LU 7\ "'l 

q—u,d 



Q 



e 2 n I d 2 k T 



P T ■ k T r j_ 2 

[xg q {x,k T ) 



M 



-D\ (z,(P T -zfc T ) 2 )]} 



(28) 



For the unpolarized TMDs ff(x,k T ), we adopt the 
results from the same spectator model calculation [4J| for 
consistency. For the unpolarized TMD FFs D\ (z,p T J, 
we assume that their px dependencies have a Gaussian 
form 



D\ (z,p 2 T ) = Dl(z) 



t(Pt) 



3 -Pt/(?t) 



(29) 



here. 



(pj.) is the Gaussian width for p T , which is cho- 
sen as (p T ) = 0.2 Gev 2 , following the fitted result in 
Ref. [67]. For the integrated FFs D\{z), we adopt the 
Krctzer parametrization [681 ] . 

To perform the numerical calculation on the asymme- 
try A s ^/ h in 7T° production at CLAS, we adopt the fol- 
lowing kinematical cuts [17] : 

0.4<z<0.7, Q 2 >lGeV 2 , W^ 2 >4GeV 2 , 

P T > 0.05 GeV, M x (eir°) > 1.5 GeV, (30) 

where M x (eir°) is the missing-mass value for the en sys- 
tem. Finally, in our calculation, we consider two differ- 
ent cases concerning the kinematical constraints on the 
intrinsic transverse momentum of the initial quarks [69j j . 
The first case is that we do not impose any constraint for 
kx in the calculation; the second case is that we consider 
the following kinematical constraints that are derived in 
Ref. H]: 



k T < (2 



X<1 



(l-2x) 2 



x)(l- 

-x) n l 



x)Q\ 



for < x < 1; 
for x < 0.5. 



(31) 



The above constraints fix the upper limit for the allowed 
kx range; thus, they can modify certain azimuthal asym- 
metries substantially, such as the the twist-3 Cahn effect 
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FIG. 3: Beam SSA A B ^ h in tt° electroproduction con- 
tributed by p x as a function of x for different Pt ranges. The 
solid and the dashed lines correspond to the results without 
and with the kinematical constraints (|31|) on fc-r, respectively. 
Data are from Ref. [lTj , the error bars for the data including 
the systematic and statistical uncertainties. 



<jf 0.10 

ui —1 
< 

0.05 



0.00 

-0.05 
0.10 

0.05 



: 0.1<xO.2 

• * • ^ ' 


0.2<x<0.3 

4 # * • 










E 0.3<x<0.4 


0.4<x<0.6 









0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 

P T (GeV) P T (GeV) 



FIG. 4: Beam SSA A a r 



in 7T° electroproduction con- 



tributed by g as a function of Pt for different x ranges. The 
solid and the dashed lines correspond to the results without 
and with the kinematical constraints (|31l) on fc-r, respectively. 
Data are from Ref. [lTj , the error bars for the data including 
the systematic and statistical uncertainties. 



analyzed in Ref. [69| . In our calculation of A s ^j h , which 
is also a twist-3 observable, we find again that the kine- 
matical constraints (|31[) modify the sizes of the asymme- 
try in certain kinematical regions. 

Figure [3] shows the results of the beam SSA A s ^f h 
for 7T° production as a function of x, compared with the 
CLAS data (full circles) measured using a 5.776 GeV elec- 
tron beam [17| . The four panels correspond to the asym- 
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metry integrated over four different Pt ranges. From the 
comparison between the theoretical curves and the data, 
one can see that our results qualitatively describe the x 
dependence of the asymmetry. Especially, for the ranges 
0.05 GeV <P T < 0.2 GeV and 0.2 GeV < P T < 0.4 GeV, 
our model calculation predicts rather flat curves, well 
agreeing with the data. We would like to point out that 
the curve for the range 0.2 GeV < Pt < 0.4 GeV is 
similar to the calculation on the g D-y term (the solid 
line in Fig. 4 of Ref. [1810 based on the models from 



Refs. 



70l | . For higher Pt ranges, deviation between 
our calculation and the data is found in the larger x re- 
gion. 

In Fig. m we display the same asymmetry from our 
calculation, but as a function of the 7r° transverse mo- 
mentum Pt- Here, the four panels in Fig. |4] correspond 
to the asymmetry integrated over four different x ranges. 
Agreement between the theoretical calculation and the 
data is found in the lower Pt region {Pt < 0.5 GeV) for 
x < 0.4. More specifically, our theoretical curves increase 
with increasing Pt in the whole Pt region within the fig- 
ure, while the experimental data increase with increasing 
Pt then approach a maximum at around Pt ~ 0.4 GeV, 
indicating that our prediction overestimates the experi- 
mental data in the region Pt > 0.5 GeV. 

From the comparison between the theoretical calcu- 
lation and the CLAS data, we conclude that the g D\ 
term can account for the beam SSA in ir° production at 
CLAS in the regions x and Pt are not large {x < 0.4 and 
Pt < 0.5 GeV). However, our model overestimates the 
data in the higher Pt and x regions. This disagreement 
might be explained by the absence of other contributions 
in Eq. (|25|) that we neglect in the calculation (such as the 
h^E), or by the possibility that the tree-level factoriza- 
tion is not suitable in this region. 

Furthermore, we also compare our calculation of the 
A ! '^/' h asymmetry for neutral pion production in SIDIS 
with the data measured by the HERMES Collabora- 
tion Q. The HERMES measurement uses a longitudi- 
nally polarized 27.6 GeV positron beam off the hydrogen 
gas target. In addition, the following kinematics are ap- 
plied in the calculation 

0.023 < x < 0.4, < y < 0.85, lGeV 2 < Q 2 < 15 GeV 2 , 
IT 2 > 4 GeV 2 , 2 GeV < P h < 15 GeV, (32) 

where Ph is the energy of the detected final-state 7r° in 
the target rest frame. In the left, central, and right pan- 
els of Fig. [5l we show the asymmetry vs z, x, and Pt 
and compare it with the HERMES data ||. It is found 
that the theoretical curves agree with the experimental 
data within the statistical uncertainty. The predicted 
z dependence of the asymmetry is rather flat, which is 
consistent with the the data in the mid-z region. The 
calculation without the constraints in Eq. (f5T|) describes 
the data in the smaller x region better than that with 
the constraints. 

Finally, we present the prediction of the asymme- 
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FIG. 5: The beam SSA for tt° production in SIDIS at 

HERMES vs z (left panel), x (central panel), and Pt (right 
panel). The solid and the dashed lines correspond to the 
results without and with the kinematical constraints (|31[) on 
kr, respectively. The thin and thick lines in the central and 
right panels correspond to the results for the ranges 0.2 < 
z < 0.5 and 0.5 < z < 0.8. Data are from Ref. []J, with 
open circles, full circles, and open squares for 0.2 < z < 0.5, 
0.5 < z < 0.8, and 0.8 < z < 1. The error bars represent the 
statistical uncertainty. 




P T (GeV) 



FIG. 6: Beam SSA of tt° at CLAS12 as a function of 

z (left panel), x (central panel), and Pt (right panel). The 
solid and the dashed lines correspond to the results without 
and with the kinematical constraints (|31|) on kr, respectively. 



try A^j h at CLAS12, which can be measured with a 
12 GeV polarized electron beam. The kinematical cuts 
for CLAS 12 applied in the calculation are [7l[ 

0.08 < x < 0.6, 0.2 < y < 0.9, 0.3 < z < 0.8, 

Q 2 > 1 GeV 2 , W 2 > 4 GeV 2 , 0.05 GeV < P T < 0.8 GeV. 

(33) 

In the left, central, and right panels of Fig. [6j we plot 
the z, x, and Pt dependencies of the asymmetry, re- 
spectively. Our calculation shows that the beam SSA 
at CLAS12 is smaller than that at CLAS, but is still 
sizable. Our result indicates that there is no obvious 
z dependence of the asymmetry in ir° production; this 
is understandable since in our approach, the same FFs 
D\(z,p T ) appear in both the numerator and the denom- 
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inator of the expression for A s ^j . Therefore, the pre- 
cision measurement of the z dependence of A s ^j h in 7r° 
electroproduction at CLAS12 can verify the role of g x in 
beam SSA. 



IV. CONCLUSION 

In this work, we have performed the calculation of the 
T-odd twist-3 TMD distribution g ± (x, k\) for the u and 
d quarks inside the proton in the spectator model with 
scalar and axial- vector diquarks. The difference between 
the isoscalar (urf-like) and isovector (ira-like) spectators 
for the axial-vector diquark is considered in the calcula- 
tion. We make use of the single-gluon exchange between 
the struck quark and the spectator to generate the T-odd 
structure. To obtain a finite result, we choose the dipolar 
form factor for the nucleon-quark-diquark coupling. We 
find that g ±u and g ±d have different x and kx dependen- 
cies. First, g ±u is positive for all x and kx regions, while 
g is negative in the small x region and turns out to be 
positive in the region x > 0.15. Second, g approaches 
to zero faster than g ±u with increasing kx- 

Using the model results, we analyze the beam SSA 
A s ^f h in semi-inclusive pion electroproduction. We ap- 
ply the Wandzura-Wilczek approximation for the FFs; 
that is, we ignore the contributions from G 1 - and E. Fur- 



thermore, we consider the specific case of ir° production, 
in which the eH^~ term should give vanishing contribu- 
tion because the favored and the unfavored Collins func- 
tions have similar sizes but opposite signs. Thus, we 
can perform the phenomenological analysis on the asym- 
metry Ail 4 " 1 in 71-0 production at CLAS and HERMES 
just based on the g ± Di term safely. The comparison be- 
tween our theoretical calculation and the data indicates 
that the g D\ term can account for the beam SSA in 7r° 
production measured by the CLAS Collaboration in the 
region x < 0.4 and Pt < 0.5 GeV, where our theoreti- 
cal curves describe the data fairly well. In addition, our 
calculated asymmetry at the HERMES kinematic region 
is consistent with the HERMES measurements after the 
error bars of the data are considered. Our study suggests 
that the T-odd twist-3 distribution g- 1 plays an impor- 
tant role in the beam SSA in SIDIS, especially in the case 
of neutral pion production. 
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